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In the effective Ising spin-1/2 antiferromagnetic chain system BaCo2V20s, the magnetic-field in- 
fluence is highly anisotropic. For magnetic fields along the easy axis c, the Neel order is strongly 
suppressed already for low fields and an incommensurate order is entered above 4 T. We present 
a detailed study of the magnetic phase diagrams for different magnetic field directions, which are 
derived from magnetization data, high-resolution thermal expansion and magnetostriction measure- 
ments as well as from the thermal conductivity. Zero-field thermal expansion data reveal that the 
magnetic transition is accompanied by an orthorhombic distortion within the ab plane. Under am- 
bient conditions the crystals are heavily twinned, but the domain orientation can be influenced 
either by applying uniaxial pressure or a magnetic field along the [100] direction. In addition, our 
data reveal a pronounced in-plane magnetic anisotropy for fields applied within the ab plane. For 
[110], the magnetic field influence on Tn is weak, whereas for magnetic fields applied along [100], 
Tn vanishes at about 10 T and the zero-field Neel order is completely suppressed as is confirmed 
by neutron diffraction data. The second-order phase transition strongly suggests a quantum crit- 
ical point being present at H ~ 10 T parallel [100], where the Neel order probably changes to a 
spin-liquid state. 



I. INTRODUCTION 

Low-dimensional quantum spin systems are of partic- 
ular interest in solid-state research due to their anoma- 
lous ground-state properties and unusual excitations, 
which are strongly influenced by enhanced quantum fluc- 
tuations. BaCo 2 V20g belongs to the material class 
AM2V2O8, which contains screw chains of MOq octa- 
hedra running along the c axis of the tetragonal crystal 
structured space group IAi/acd with lattice constants 
a = 12.444(1) A and c = 8.415(3) A. In the ab planes, 
the chains are separated by nonmagnetic A 2+ (A = Sr, 
Ba) and V 5+ ions. Depending on the transition metal 
ion M 2+ (M = Cu, Ni, Co, Mn), the various mem- 
bers of y4M 2 V 20g represent spin-chain materials with 
different spin quantum numbers and anisotropics^^. Ac- 
cording to Hund's first rule, the Co 2+ ions realize the 
5 = 3/2 high-spin state and the lower-lying ti g or- 
bitals of the crystal-field split 3d states are only par- 
tially filled, resulting in a finite effective orbital moment 
/ = 1 for regular CoOg octahedra. In BaCo2V20s, the 
CoOg are compressed along the c axis, which causes a 
pronounced Ising anisotropy of the total magnetic mo- 
ment with c as the easy axis. This Ising anisotropy is 
experimentally confirmed by low-temperature high-field 
magnetization datcP which reveal a strong anisotropy of 
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both, the saturation field H^ at sa 0.5 H^. and the sat- 
uration magnetization M^ t sa 1.5 as well as by 
temperature-dependent magnetic-susceptibility measure- 
mental in the intermediate temperature range yielding 
x" c (T) sa 2 X ±C (T) for T > 5.5 K. Moreover, x" c (T) 
shows a pronounced maximum around 35 K, which sig- 
nals short-range magnetic correlations reflecting the one- 
dimensional nature of the magnetic subsystem due to 



the dominant magnetic exchange along the CoOe screw 
chains. Thus, BaCo2V20g represents a model system of a 
(quasi-)one-dimensional antiferromagnetic effective Ising 
5 = 1/2 chain that is described by the XXZ hamiltonian 

ft = E J { s t s Ui + e ( s i s ?+i + s i s !+i) } + gvBSiH . 
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From an analysis of the low-temperature magnetization 
data, J ~ 65 K and e ~ 0.46 as well as the anisotropic g 
factors <?" c ~ 6.3 and g ±c ~ 3.2 have been derivecP. 

Due to finite inter-chain couplings, the system shows 
long-range antiferromagnetic order below Tn ~ 5.5 K 
with a propagation vector &af = (1,0,0) and the spin 
direction along the ±c axis. Nearest-neighbor (NN) spins 
along c are oriented antiparallel to each other, while 
within the (001) planes, NN spins point either in the same 
or in the opposite directiorpSHEl As a consequence, the 
magnetic symmetry is less than tetragonal and the or- 
dered phase typically contains different domains, which 
are rotated by 90° and translated by c/4 with respect to 
each other^. The magnetic-field vs. temperature phase 
diagram of BaCo 2 V 2 8 has been studied by several tech- 
n jq ue ^Mlll an d a pronounced anisotropy has been re- 
portedl^LLSI For H\\c, Tn rapidly decreases in the field 
range up to about 4 T and above this field an incommen- 
surate (IC) magnetic phase with a weakly field-dependent 
transition temperature Tic ~ 1 K develops. Moreover, 
experimental evidence for additional phase tra nsitio ns in 
the region of the IC phase has been reported 17 * 18 ^. For 
H _L c, a rather weak decrease OT^/dH i=a —0.1 K/T has 
been reported up to a maximum field of 9 T in Ref. 1131 
whereas a much stronger decrease of dT^/dH sa —1 K/T 
for 8 < H < 10 T has been found in Ref. El However, 
in most of the published studies for H _L c, the explicit 
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direction of t he magneti c field within the ab plane has 
not been giver^M 1 3 1 17 1 20 " ^. 

In this report, we present a detailed study of the 
magnetic phase diagrams of BaCo2V20s for differ- 
ent directions of the magnetic field. Based on high- 
resolution thermal-expansion measurements by capaci- 
tance dilatometry, we find that the magnetic ordering 
in zero magnetic field is accompanied by small struc- 
tural distortions, which break the tetragonal symmetry. 
The distortions arise from the magnetoelastic couplings 
within the ab plane and the crystals are twinned accord- 
ing to the magnetic domain structure. For H\\c, our data 
confirm the existing phase diagram in the low-field range 
and reveal a nonmonotonic field dependence of Tic to- 
wards higher fields. For H _L c, we find an additional 
anisotropy, which has been overlooked so far and resolves 
the apparent discrepancy between the phase diagrams of 
Refs. [13] and [THl For a magnetic field applied along the 
[110] direction, we observe, in fact, the above-mentioned 
weak decrease of dT^/dH ss —0.1 K/T up to our max- 
imum field of 16 T, but for ff||[100] (respectively || [010] 
of the twinned crystals), a much stronger decrease of Tn 
occurs (in agreement with the data of Ref. H9|) and our 
data suggest a quantum phase transition at H c ~ 10 T. 
This additional in-plane anisotropy of the magnetic-field 
influence has been also reported by a very recent mag- 
netization study'™ Moreover, our measurements reveal 
that the crystal can be (at least partially) detwinned by 
applying either uniaxial pressure or a magnetic field along 
[100]. The paper is organized as follows: Details of the 
single crystal growth and the experimental setups are de- 
scribed in section [TTJ while the experimental results are 
presented and discussed in section [TTTJ which is split into 
4 subsections discussing the zero-field measurements and 
the results for the different field directions -ff ||c, -ff |j [110], 
and -ff||a, respectively. 



II. EXPERIMENTAL 

Single crystals of BaCo2V20s have been synthesized 
by spontaneous nucleation as well as by the floating- 
zone method. In both cases, BaCo 2 V 2 08 was synthe- 
sized by a solid-state reaction using a mixture of BaCC>3 
(99+ % Merck), Co 3 4 (99.5 % Alfa Aesar) and V 2 5 
(99.5 % Strem Chemicals) in a 1:|:1 molar ratio. A 
four-mirror image furnace (FZ-T-10000-H-VI-VP, Crys- 
tal Systems Inc.) was used for the single-crystal growth 
by means of the floating-zone technique (in analogy to 
Ref. I22p . Atmospheric air at ambient pressure served as 
reaction atmosphere. Fig. [T] shows photographs of both 
crystals. The single crystal grown by spontaneous nucle- 
ation is of rectangular shape with dimensions of about 
1.2 x 1.3 x 2.1 mm 3 . The long axis is parallel to the c 
direction and the flat surfaces are (110) planes. 

The magnetization has been measured in a vibrat- 
ing sample magnetometer (PPMS, Quantum Design) for 
temperatures above 2 K up to a maximum magnetic 




FIG. 1. Single crystals of BaCc^V^Og obtained by the 
floating-zone method (b) and by spontaneous nucleation (c). 
The arrow in (c) marks the magnetic easy c axis and the flat 
surface is a (110) plane. Panel (a) summarizes basic features 
of the crystal and magnetic structure. Right- and left-handed 
screw chains of CoOe octahedra are running along c. For 
clarity, four of these chains are also shown without the O 2- 
ions and the dominating NN exchange J along the chains, see 
Eq. |l]), is depicted by solid lines. The dashed and dotted 
lines connect NN Co 2+ ions in the ab planes at z — 1/8 and 
z = 3/8, respectively, and represent NN inter-chain couplings 
J^ N , which compete with the corresponding NNN inter-chain 
couplings Jjvjviv- Panel (d) displays the geometrical arrange- 
ment of the Co 2+ ions in the ab planes. Due to the competing 
inter-chain couplings Jjvjv an d Jnnn^ this arrangement rep- 



resents a mo re comp lex version of the Ji 
square lattice 



J 2 model on a 



field of 14 T. Temperature or magnetic-field induced rel- 
ative length changes AL(T, H)/L, i.e. thermal expansion 
and magnetostriction, have been determined by a home- 
built capacitance dilatometer in the temperature range 
from about 250 mK to 10 K in magnetic fields up to 
16 T. Moreover, zero-field thermal-expansion measure- 
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ments have been performed in another dilatometer, which 
allows to easily vary the weak uniaxial pressure applied 
by the springs connecting the moveable and fixed ca- 
pacitor plates. The corresponding expansion coefficients 
a(T) = l/LdAL(T)/dT and \(H) = l/LdAL(H)/dH 
have been obtained numerically. In finite fields, the 
length change has been determined parallel to the direc- 
tion of the applied field, i.e. L,-||_ffj. The thermal conduc- 
tivity has been measured by the standard steady-state 
method from ~ 0.3 to f5 K. The temperature gradient 
was produced by a heater at one end of the sample and 
measured by two matched Ru02 thermometers. The di- 
rections of the heat current were either along the c axis or 
the [f f 0] direction of the crystal of Fig. [TJ; and magnetic 
fields up to ~ 15 T have been applied along c, along a, 
and along the [110] direction. Neutron diffraction studies 
were performed at T ~ 50 mK in magnetic field up to 
12 T on the CEA-CRG single-crystal D23 diffractome- 
ter, operated in the normal beam mode with an incident 
wavelength of 1.2765 . The sample of about 250 mm 3 
was cut from the crystal rod obtained by the floating- 
zone method and was aligned with the a axis vertical. 
It was mounted in a 12 T vertical cryomagnet equipped 
with a dilution insert, with the field applied along the o 
axis. 



III. RESULTS AND DISCUSSION 
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FIG. 2. (color online) Zero-field thermal expansion of the 
tetragonal at axis of BaCo2V2 08- The different data sets 
in (a) were obtained on two different setups and samples, 
while those in (b) were measured during the same run, but 
with varying uniaxial pressure Pi applied along at . All length 
changes AL — L(T) — Lq are related to Lq = L(7 K,Pi). The 
process of a partial detwinning by uniaxial pressure or in a 
magnetic field is schematically sketched in the lower panel. 



A. Zero-field thermal expansion 

Fig.[2]displays thermal-expansion data along the a axis 
of BaCo2V208 measured in zero magnetic field. The left 
panel shows AL(T)/L obtained on two different samples 
in different setups. In both cases, the magnetic ordering 
is accompanied by spontaneous strains, but these strains 
are of different signs. Because tetragonal symmetry re- 
quires an isotropic thermal expansion within the ab plane, 
this observation is an unambiguous proof that the tetrag- 
onal symmetry is broken in the Neel-ordered phase. This 
can be naturally explained by assuming that the antifer- 
romagnctic ordering causes a small orthorhombic split- 
ting a ^ b below T N and that the crystals are twinned 
with respect to the orientations of a and b. Consequently, 
the measurements of Fig. [2^, reflect different superpo- 
sitions of the a- and &-axis expansions due to different 
twin ratios in both runs. In order to check this, we per- 
formed a series of thermal-expansion measurements along 
the tetragonal at axis of one crystal for varying uniaxial 
pressure along this axis. The uniaxial pressure p t arises 
from the spring which fixes the sample and the moveable 
plate to the fixed capacitor plate and increases with in- 
creasing base capacitance Co <x 1 /do , where do is the base 
distance between both capacitor plates. As the spring of 
the actual setup is rather soft, the resulting pi for a sam- 
ple with a cross section in the mm 2 range can be varied 
around ps 0.5 MPa. 

As is shown in Fig. [2b, increasing pi from 0.28 MPa 



to 0.76 MPa is indeed sufficient to switch from a sponta- 
neous expansion below Tn to a spontaneous contraction. 
In this low-pressure range, the magnitude of the sponta- 
neous strain continuously varies with increasing pi, what 
is naturally interpreted as follows. Below Tn, the tetrago- 
nal at axes spontaneously develop an orthorhombic split- 
ting, e.g. a Q > a t > b a , but due to a twinning of the 
ferroelastic domains, the measured macroscopic length 
change is an average of a (T) and b (T). With increas- 
ing pi along one of the tetragonal a t axes, the amount of 
& Q -oriented domains along this direction is continuously 
enhanced until a single-domain crystal would be reached. 
This process is schematically sketched in Fig. [2b, but 
the uniaxial pressure applied here is not sufficient to 
reach this single-domain limit. The detwinning process 
should be directly visible as a function of Pi below Tn, 
but similar to the low-temperature magnetization [polar- 
ization] loops M(H)[P(E)] in ferromagnetic [ferroelec- 
tric] materials, pronounced hysteresis effects may occur. 
Therefore (and because of technical reasons in our setup) , 
we performed "pressure-cooled" measurements meaning 
that the uniaxial pressure has been changed well above 
Tn at ~ 10 K, then the sample has been cooled to ~ 3 K 
and the data of Fig. [2}d have been collected during the 
subs equent heating run. As will be discussed below (Sec- 
tion III C 2 ) , an analogous detwinning is also achieved by 



applying a magnetic field along a t and in this case, we 
can also observe the typical hysteresis effects of zero-field 
cooled measurements. 



The structural symmetry reduction in the N eel-ordered 
phase of BaCo2V2 08 has not been discussed in any of the 
reported diffraction studies^ 12 * 16 1 , what is most proba- 
bly due to the fact that the orthorhombic splitting is 
very small. As shown in Fig. |2j the averaged macro- 
scopic strains AL/L are of the order of 10 -6 and as will 
be discussed below, we expect that even in a completely 
detwinned sample AL/L does not significantly exceed 
10~ 4 . Thus, it will be hard to detect such small changes 
by diffraction techniques. Nevertheless, our finding is 
supported by the m agnetic structure obtained by neu- 
tron scatterin g* 11 ! 12 !. As is shown in Fig. [l] NN spins in 
the ab plane are oriented parallel to each other along 
one tetragonal at axis, but antiparallel to each other 
along the other a t axis. This magnetic structure typi- 
cally arises from magnetic frustration, e.g. a competition 
between antiferromagnetic NN interactions Jnn along 
at and antiferromagnetic interactions Jnnn between the 
next-nearest neighbors along the diagonal (110) direc- 
tions^ see Fig. [l] and the corresponding discussions in 
Refs. HH I25H291 Here, the notations are used be- 
cause we are discussing inter-chain interactions, whereas 
the dominant magnetic intra-chain interaction J of the 
hamiltonian ([I]) acts along the chain direction. The com- 
petition between J^ N and Jnnn results in a perfect 
magnetic frustration for a tetragonal symmetry where 
both at directions are equivalent, but it can easily be 
lifted by magnetostriction if there is a finite magnetoe- 
lastic coupling, i.e. if J^ N and/or J^nn change with 
lattice distortions. One of the simplest scenarios would 
be that Jnn decreases with increasing NN distance. In 
this case, a small orthorhombic distortion a Q > at > b Q 
lifts the magnetic frustration and couples the magnetic 
and lattice structures such that NN spins are parallel to 
each other along a and antiparallel along b Q . Note that 
this oversimplified model could already qualitatively ex- 
plain the above observation, but at present, neither the 
sign nor the magnitude of the change of Jnn as a func- 
tion of the distance is known and, moreover, Jnnn can 
also depend on the lattice distortions. 



B. Magnetic field parallel to the c axis 

The thermal-expansion and magnetostriction data 
along the chain direction c for different magnetic fields 
H\\c, i.e. along the Ising axis, are summarized in Fig. [3] 
The Neel ordering in low fields is accompanied by a 
pronounced spontaneous contraction (Fig. [3^,) and the 
field-dependent Tn(H) can be easily monitored by con- 
sidering the shift of the corresponding a(T) anomalies 
(Fig. [3}d). The shape and magnitude of these anomalies 
hardly change up to about 3 T, but with further increas- 
ing field the anomalies rapidly drop in size and drastically 
broaden. As shown in Fig. [3j;, the position as well as the 
magnitude of these anomalies non-monotonically vary 
with increasing H. The drastic change of the anoma- 
lies around 4 T is related to a transition from the Neel- 
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FIG. 3. (color online) Thermal expansion (left) and magne- 
tostriction (right) of the c axis of BaCo2V20s for magnetic 
fields H\\c. In the upper panels the relative length changes 
(related to L = L(7 K, H) and L = L(T, T), respectively) 
are shown, whereas the lower panels display the derivatives of 
AL/L with respect to T or H. Panel (c) is an enlarged view 
of the low-temperature range of (b); the different curves in 
(e) are measured with increasing (•) or decreasing (o) field. 



ordered state to an incommensurate high-field phase as 
has been established by diffraction studief^^]. The cor- 
responding length change have also been studied by low- 
temperature magnetostriction measurements AL{H)/ L, 
see Fig. [3}i. Below about 1.5 K, the length changes at 
« 4 T are almost discontinuous, and depending on the 
field-sweep direction, there is a clear hysteresis confirm- 
ing the proposed l st -order character of this transition 
(Fig. [3j3) . Above 1.5 K, the magnetostriction curves are 
more continuous and the hysteresis vanishes, i.e. the 
transition is of 2 nd order. 

At least qualitatively, all the features of the field- 
and temperature-dependent length changes AL(H,T) of 
BaCo2V208 strongly resemble the spontaneous strains 
observed in the spin-Peierls compound CuGeOj^nHlH. 
In both materials, the spontaneous strains are large in 
the respective low-field commensurate phase and become 
considerably reduced when an incommensurate high-field 
phase is entered. Moreover, the temperature-induced 
transitions into the incommensurate phases are broad- 
ened compared to the transitions into the commensurate 
phases. For doped CuGeC>3 , it has been observed that the 
degree of broadening is enhanced with increasing doping 
concentration suggesting that the incommensurate phase 
is more severely disturbed by disorder than the commen- 
surate phased Because the field-induced broadening in 
pure BaCo2V2 0g is already more pronounced than that 
in weakly doped CuGeOa, one may speculate that the 
rather complex structure of BaC^A^Os could be very 
sensitive to defects and this could be the reason for the 
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FIG. 4. (color online) Phase diagram of BaCo2V208 for 
H\\c. While the transitions to the Neel-ordered phase are 
very sharp, the vertical bars in the higher-field range mark 
broad crossover regions to the incommensurate phase. 



drastic broadening of the transitions to the incommensu- 
rate phase. However, the spin-Peierls order in CuGe03 
can obviously not be directly compared to the Neel order 
of BaCo2V2 0g where the microscopic nature of the high- 
field phase is not yet clear (see below). Nevertheless, the 
mere fact that there is a pronounced spontaneous con- 
traction of c in the Neel-ordered phase establishes a sig- 
nificant magnetoelastic coupling. Again, we consider the 
most simple scenario of a NN interaction J, see Eq. ([lj, 
which increases with decreasing NN distance. In this 
case, the spontaneous contraction Ac below Tn would 



result from a gain in magnetic energy (£"„ 



-Ac 



which ovcrcompensates the corresponding loss of elastic 
energy (E e i oc Ac 2 ). Together with the above discussion 
of the orthorhombic distortions a and b , we can, thus, 
summarize that the perhaps most simple model of anti- 
ferromagnetic NN intra-chain (J) and inter-chain ( J^ N ) 
interactions can explain the observed lattice distortions 
below Ttf. We emphasize, however, that this scenario is 
oversimplified because (i) the considered NN couplings 
J l \\i do not only depend on the corresponding longitu- 
dinal distortions but also on the distortions in the 
transverse directions j, k _L i and (ii) additional, e.g. 
NNN, couplings and their dependencies on distortions 
along k may also be important for a quantitative de- 
scription. 

In Fig. [4j we present the corresponding phase diagram 
of BaCo2V2C>8 for magnetic fields along the c axis. The 
phase boundaries are obtained from the positions of the 
extrema of a(T) and X(H). Additional data points based 
on our magnetization and thermal-conductivity measure- 
ments are also included. In all cases, the anomalies at the 
transitions to the Neel-ordered phase are rather sharp, 
while those related to the incommensurate phase are 
very broad. Because the high-field anomalies rather in- 
dicate a broad crossover behavior instead of well-defined 



phase transitions, we marked the positions of the max- 
ima of a(T) in the phase diagram of Fig. [1] as vertical 
bars representing the anomaly width. The phase dia- 
gram agrees well to those presented previously by us- 
ing other techniques' 12 ' 14 ' 15 ' 17 ' 18 ' ^ . In the low-field 
range of the Neel-ordered phase all these reports essen- 
tially agree with each other, but some discrepancies are 
present in the high-field range. For example, additional 
features have been observed in sound-velocity measure- 
ments in the field range close to the triple pointP^, but 
corresponding effects are neither seen in the other litera- 
ture data nor in our measurements. Most of the previous 
reports could follow t he bo undary to the high-field phase 
only up to ~ 10 iCnHHHltni anc i on i y reC ently, the field 
range has been extended up to 14 T by an NMR studjfiSl. 
These measurements revealed a minimum of the phase 
boundary around 8.6 T and suggested that this mini- 
mum is related to a rearrangement concerning the mag- 
netic inter-chain ordering. Theoretical studies on the ba- 
sis of bosonization combined with mean-field theory for 
the inter-chain interaction suggested an alternative pos- 
sible origin for a non-monotonic phase boundary, namely 
a possible change from the incommensurate phase with 
spins parallel to c towards a transverse stagg ered p hase 
with spins oriented (mainly) perpendicular to <j 33 ' 34 l De- 
spite the drastic broadening of the high-field anomalies, 
our data clearly support the occurrence of a minimum of 
the phase boundary around «8T, but our macroscopic 
data can neither confirm nor refute any of the proposed 
microscopic changes of the magnetic structure. In our 
low-temperature magnetostriction data (Fig. |3j;) , we do 
not observe any indication for the occurrence of a field- 
induced phase transition around 8 T. However, in view 
of the very small and broad anomalies of a(T) at the 
boundary to the incommensurate phase, one might also 
expect only extremely shallow anomalies at a possible 
rearrangement of the magnetic structure. 



C. Magnetic field perpendicular to the c axis 

Up to now, the field direction perpendicular to the easy 
c axis of BaCo2V2C>8 seemed to be of comparatively lit- 
tle interest. As described above, low-temperature M(H) 
curves^ yield a saturation field of about 40 T and the re- 
ported phase boundaries in the region of 9 T are contra- 
dictorjUSEq, In both studies, the explicit field direction 
within the ab plane has not been specified, but our mea- 
surements reveal a significant anisotropy between if || at 
and || [110]. Such an in-plane anisotropy is in accor- 
dance with the tetragonal symmetry of the paramagnetic 
high-temperature phase and as shown above, the sym- 
metry in the Neel-ordered phase is even reduced to or- 
thorhombic; although for macroscopic probes this is (at 
least partly) masked by the twinning of the crystals. 

Fig. [5] compares representative magnetization data of 
BaCo 2 V 2 8 for H\\a t and #||[110]. For H = 1 T, the 
M (T) curves for both field directions only differ by about 
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FIG. 5. (color online) Left: T-dependent magnetization of 
BaCo2V2 08 for different magnetic fields applied either along 
[110] (top) or along the tetragonal at axis (bottom). Right: 
Magnetization curves M(H) (top) and the corresponding 
dM/dH (bottom) for both field directions, measured at 2.3 K. 



10 % and show a sharp kink around T/v — 5.5 K. With 
increasing field, however, a pronounced anisotropy de- 
velops: Increasing the field to 10 T for if||[110] causes 
a continuous suppression of M/H and a weak decrease 
of Tn, whereas for H\\a t , there is little change of M/H 
above Tn , but Tn itself is so strongly suppressed that it 
is no longer visible in the 10 T curve. The right panels of 
Fig. [5] compare the M(H) curves measured at T = 2.3 K 
for both field directions. While M{H) for if || [110] just 
continuously increases (with a slight downward curva- 
ture) up to the maximum field of 14 T, the M(H) curve 
for if 1 1 at has a clear anomaly at ~ 9.2 T signalling a 
phase transition as is best seen in the corresponding dif- 
ferential susceptibility dM/dH. A very recent high-field 
magnetization study has revealed that this phase transi- 
tion is not related to a saturation of M and has proposed 
that its microscopic origin lies in the structural arrange- 
ment of the CoC>6 screw chainiPSl (see below). 

In the following, the phase diagrams obtained for the 
two in-plane field directions are discussed separately. 
Again, the corresponding phase boundaries are derived 
by measurements of the thermal expansion and magne- 
tostriction, which do not only have a very high resolution 
but also yield additional information about the magne- 
toelastic coupling and a field-induced reorientation of the 
domains. 



Magnetic field parallel to the [110] direction 



Representative measurements of AL[ 110 ] (T)/Lni ] for 
different fields TT 1 1 [1 10] are displayed in Fig. [6j In 
zero field, a spontaneous contraction develops below Tm, 
which is weakly enhanced on increasing the magnetic 
field, while Tn, derived from the maxima of the corre- 
sponding a[ 110 ](T) (not shown), continuously decreases 
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FIG. 6. (color online) Thermal expansion AL(T)/L measured 
along [110] for various magnetic fields i/||[110] (left) and the 
corresponding H — T phase diagram (right). 



with a slope dT N /8H ~ -0.1 K/T. Thus, our phase 
boundary for if || [110] (Fig. [6]) very well agrees to the 
H ±c - T phase dia gram based on specific-heat dataP^ 
and we suppose that those data were measured with 
if -1 - !! [110]. Such an orientation appears natural because 
the crystal of Ref. [T31 was obtained by spontaneous nu- 
cleation and our crystal grown by this method has flat 
(110) surfaces; see Fig. [TJ;. 

As discussed, the Neel ordering causes an orthorhombic 
splitting a > a t > b and the crystal is twinned below 
Tn- Such a twinning does, however, not influence the ex- 
pansion along [110] because both types of domains are ro- 
tated with respect to each other by 90° around the c axis 
and the relative length change along the diagonals of each 
domain is given by Aiiio/iiio = (Ao /a + Ab /b ) /2. 
The individual Aa /a and Ab /b are of opposite signs 
(Fig. [2]) , but we expect their absolute values to be of sim- 
ilar magnitude and not too much larger than the mea- 
sured AL110/T110 ~ 1.5 • 10~ 5 . The corresponding or- 
thorhombic splitting (a a — b a ) / (a D + b D ) should thus be 
in the 10 -4 range making it difficult to be detected by 
standard diffraction techniques as already stated above. 



2. Magnetic field parallel to the a axis 

Now we turn to the configuration with H applied 
along one of the tetragonal a t axes. The thermal- 
expansion data measured along this direction are pre- 
sented in Fig. [7j In this setup, the domain structure 
below Tn is such that a small spontaneous contraction is 
present in zero field. On increasing the magnetic field, 
the magnitude of the spontaneous contraction rapidly de- 
creases and around 3 T, it even changes sign and turns 
into a spontaneous expansion, which is of maximum size 
around 6 T and finally completely vanishes above 10 T. 
The field dependence of Tn is best monitored by con- 
sidering the shift of the maxima/minima of the corre- 
sponding thermal-expansion coefficient a(T) shown in 
Fig. [3]d. At first glance, the sign change around 3 T 
might be interpreted as an indication for different mag- 
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FIG. 7. (color online) Thermal expansion (left) and magne- 
tostriction (right) of the tetragonal at axis of BaCo2V208 for 
magnetic fields H\\at- The upper and lower panels display 
AL/L and the corresponding derivatives with respect to T 
and H, respectively. For clarity, the different curves in pan- 
els (b) and (d) are shifted with respect to each other. Around 
3 T there is a sign change of AL(T) due to a field-induced par- 
tial detwinning and for the same reason the field-dependent 
data depend on the field-sweep direction (marked by arrows) 
in the low-field range. 



netic structures above and below that field. However, 
we already know that the crystal is twinned and a more 
natural explanation arises from a magnetic-field induced 
(partial) detwinning by assuming that the magnetic field 
favours those domains which are oriented with the longer 
a a axis parallel to H. Thus, if applied along the same 
a t axis, a magnetic field and uniaxial pressure would op- 
pose each other because H favors the "long" domains be- 
ing aligned along -ff||a f whereas uniaxial pressure prefers 
the "short" domains along its direction. However, in a 
transverse configuration, i.e. when applied along two 
orthogonal a t axes, magnetic field and pressure would 
support each other because both favor the same type of 
domain orientation; see Fig. [2j As explained above, the 
AL(T,pi)/L data for uniaxial pressure (Fig.J^) were ob- 
tained in a "pressure-cooled" mode and, analogously, the 
AL(T, H)/L curves of Fig. [7^, were measured in the "field- 
cooled" mode, i.e. H has been varied only above Tn- 

The right panels of Fig. [7] display some AL(H)/L 
curves measured at constant T after the sample has been 
cooled in zero field. At T = 0.31 K, AL(H)/L shows two 
anomalies on increasing the field: a step-like increase at 
~ 6 T and a slope change at ~ 10 T. The slope change is 
reversible, whereas the step-like change is absent in the 
field-decreasing run. At a higher T = 3 K, broadened 
but still hysteretic anomalies occur in AL(H)/L around 
2 T for both field-sweep directions, while at ~ 8.5 T, 
the analogous reversible slope change is observed as in 



the low-temperature curves. The slope changes at higher 
fields correspond to a field-driven crossing of the 2 nd - 
order phase boundary, which is already known from the 
temperature-dependent measurements of the left panels 
of Fig. [7] The anomalies at lower fields show a similar 
hysteresis behavior as it is known from hysteresis loops 
of ferromagnets and the more or less step-like increase 
on increasing field is, thus, naturally explained by the 
field-induced domain orientation as discussed above. 

Note, however, that in contrast to a ferromagnet, the 
field-induced orientation of the domains in BaCo2V20s 
cannot be reversed by reversing the magnetic-field di- 
rection. Instead one would have to apply either a mag- 
netic field along the perpendicular at direction or, al- 
ternatively, one could apply uniaxial pressure along the 
direction of the original field. The possibility of manipu- 
lating a certain type of "ferro-" domains, e.g. ferromag- 
netic, ferroelectric and/or ferroelastic, by different types 
of external fields, i.e. magnetic or electric fields or uniax- 
ial pressure/strain, is of strong potential interest in the 
context of data storage, and during the last decade it 
has attracted a lot o f atte ntion in the field of so-called 
multiferroic material d 35 * 36 l Typical multiferroic ordering 
phenomena arise from coupled (anti-)ferromagnetic and 
ferroelectric order parameters and, in most cases, a more 
or less complex magnetic structure yields the primary 
order parameter, while the ferroelectric polarization is a 
secondary (or improper) order parameter. The coupling 
of both order parameters, then, allows the manipulation 
of the magnetization by an electric field or, vice versa, 
changing the polarization by a magnetic field, and these 
changes are most pronounced if there are switchable do- 
mains. In the present case of BaCo2V2 0s, the primary 
order parameter stems from the (antiferro-)magnetic or- 
der while the spontaneous strain is a secondary order pa- 
rameter and, due to the finite order-parameter coupling, 
the corresponding domains can be reoriented either by a 
magnetic field or by uniaxial pressure. 

Fig. [8] displays representative thermal-conductivitiy 
data of BaCo2V2C>8 for different magnetic field direc- 
tions. In all these measurements, the heat current has 
been applied along c, i.e. parallel to the spin-chain di- 
rection. One aim of these measurements was to study 
whether there is a sizable magnetic contribution to the 
total heat transport as it is foun d in s everal Heisenberg 
spin-l/2-chain or ladder material d 37 * 38 l A standard way 
to separate a quasi-lD magnetic contribution from the 
phonon heat transport is to consider the difference of the 
thermal conductivity along (k") and perpendicular (fir 1- ) 
to the chain direction. However, our low-temperature 
data of k (not shown) obtained with a heat current 
along [110] are almost identical to the (T) of Fig. [8] 
The lack of anisotropy between k" and is a strong in- 
dication against the presence of a sizable magnetic heat 
transport in BaCo2V2 0s. This finding has recently been 
published already in Ref. [19] by Zhao et al. who have 
drawn exactly the same conclusion. As is seen in Fig. [HJ 
the zero-field k(T) has a sharp minimum at Tjy — 5.5 K, 
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c of BaCo2V208 for various magnetic fields applied along 
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Kc C (H) I kI c (0) for 7? || at at various temperatures. 



which signals strong phonon scattering by magnetic fluc- 
tuations around the ordering transition and following the 
minimum of k for different magnetic fields (and field di- 
rections) reveals the respective phase boundaries. The 
k(T) data for different field directions fully confirm the 
strongly anisotropic phase diagrams of BaCo2V20g for 
H\\a, if || [110], and H\\c as discussed above. Moreover, 
our «(T, if) data for if||a and if||c (Figs. [8^, b, and d) 
very well agree to the respective measurements of ft(ff , T) 
for if _L c and if || c of Ref . [TSJ Because the additional in- 
plane anisotropy was unknown, Zhao et al. interpreted 
the deviation of their if - Lc -phase boundary from the pre- 
viously known 7J- Lc -phase boundary^ as indication for 
the occurrence of a new phase around 10 T. This dis- 
crepancy is, however, due to a mixing up of the different 
phase diagrams for if | [110] and if ||a and is fully resolved 
by considering the very different results of k(T, if) for the 
different in-plane field directions if _L c, see Figs. [8^, and 
[8b. 

The phase diagram of BaC^A^Og for magnetic fields 
if || a is given in Fig. [9} The shape of this phase boundary 
is consistently confirmed by all the different probes. The 
small quantitative deviation of the critical fields obtained 
by one set of thermal-conductivity measurements com- 
pared to the other data most probably arises from a slight 
misalignment of the magnetic-field direction. We men- 
tion that a similar strong suppression of Tn for if _L c has 
been reported for the closely related material SrCo2V2 0g 
in Ref. [351 Because SrCo 2 V20g and BaC^A^Og have 
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FIG. 9. (color online) Phase diagram of BaCo2V2 08 for ff||at. 
The dashed line marks the region of the field-induced reori- 
entation of magnetic domains. 



essentially the same structure and also almost the same 
Neel temperature in zero field, their magnetic subsys- 
tems should be very similar and we expect that a study 
of the in-plane anisotropy of SrCo2V20g would reveal 
very similar results as obtained here for BaCo2V20g. As 
indicated by the dashed line in Fig. [9j the partial field- 
induced detwinning of the sample causes some hysteresis 
behavior, which not only depends on the actual path in 
the if — T plane but is also sample dependent because it 
is probably influenced by internal strains and/or defects. 
Concerning the actual phase boundary, however, we did 
not observe hysteresis effects in any of our measurements 
identifying the corresponding transition as a continuous 
2 nd -order phase transition. Thus, our data strongly sug- 
gest that BaCo2V20g has a quantum critical point at 
H& t * 10 T. 

As can be seen in Fig. [7]d, the thermal expansion co- 
efficient a also changes sign at 10 T, which is a typ- 
ical signature of a quantum phase transitio n 140 * 41 *. For 
purely one-dimensional magnets, this sign change can 
even take place in form of a pole, i.e. a(T — ¥ 0) di- 
ver g e JHH3] whereas for higher-dimensional systems, such 
a sign change via a divergence is only present for the 
Gruneisen ratio a/C (C is th e spec ific heat) and the 
magnetocaloric effect d In T j e>j j 41 l 44 l EI] In the present 
case of BaCo2V20g, we are most probably dealing with 
a three-dimensional quantum phase transition from the 
Neel-ordered phase to a spin-liquid state. Very recently, 
this transition has also been observed by high-field mag- 
netization data and the strong in-plane magnetic-field 
anisotropy of BaCo2V2 0g has been explained by a struc- 
tural peculiarity of the Co 2+ screw chains arising from a 
slight distortion of the CoC>6 octahedral. According to 
this model, a homogenous magnetic field if || a is accom- 
panied by a staggered transverse field component via a 
staggered g tensor, whereas no such transverse compo- 
nent occurs for if || [110]. As a consequence, this model 
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FIG. 10. (color online) Intensities of several magnetic Bragg 
peaks of the Neel-ordered phase of BaCo2V20s as a function 
of magnetic fields H\\a t at T ~ 50 mK. The peaks of panel 
(a) arise from one typ e of domains while those of (b) repre- 
sent the other type^H The different intensities obtained with 
increasing or decreasing magnetic field, shown as solid and 
open symbols, respectively, proof a field-induced reorienta- 
tion of the magnetic domains in the ordered phase. 



predicts a collapse of the Neel order at already rather 
small fields of ~ 7 T applied along the a direction. 

As is shown in Fig. [TUJ this hypothesis is verified by 
our low-temperature (T ~ 50 mK) single-crystal neutron 
diffraction data. The field dependent intensities of sev- 
eral magnetic superstructure reflections of the zero-field 
Neel-ordered phase collapse at H^ it ~ 10 T. However, 
the suppression of the magnetic intensity is not homoge- 
neous and there is a significant hysteresis. Moreover, all 
reflections show a clear anomaly around 4 T, whose origin 
remains unclear. Based on the present data, it cannot be 
excluded that this anomaly is related with the domain 
structure and extinction effects. As discussed above, the 
zero-field Neel phase contains two types of magnetic do- 
mains, which are represented by the different magnetic 
superstructure peaks shown in Fig. |10fr, and |T0p , respec- 
tively (see Ref. [T^)) . On increasing the magnetic field, 
the the peak intensities of panel (a) rapidly decrease al- 
ready above ~ 9 T, while those of panel (b) only drop 
close to 10 T. Such a different field dependence is nat- 
urally explained by a field-induced domain reorientation 

close to the critical field H^ it ~ 10 T. Moreover, the fact 
that one type of domains is favored by the magnetic field 
H\\a t also explains why the field-decreasing run results in 
a different intensity distribution as compared to the orig- 
inal zero-field data. The suppression of the intensities of 
all measured peaks above H^ it ~ 10 T implies that ei- 
ther the high-field phase is indeed a spin-liquid state or 
that another type of magnetic order settles above H^ it . 



Because the present data set does not yet allow an un- 
ambiguous conclusion, an additional neutron scattering 
study of the high-field phase is needed to look whether 
some magnetic signal appears elsewhere in the reciprocal 
lattice or not. 

IV. SUMMARY 

We have presented a detailed study of the magnetic 
phase diagrams of BaCo2V20g for field directions par- 
allel and perpendicular to the easy axis c. For a mag- 
netic field along the easy axis, our data essentially con- 
firm previous results, in particular the non-monotonic 
field dependence of the transition temperature to the in- 
commensurate phased The phase diagram for magnetic 
fields applied within the ab plane is more complex than 
originally expected. The zero-field Neel ordering is only 
weakly suppressed by a magnetic field along [110]. How- 
ever, applying the magnetic field along the tetragonal 
a t axis, the phase diagram shows a strong field depen- 
dence and the Neel temperature vanishes above a criti- 
cal magnetic field H'^ it ~ 10 T. Our thermal-expansion 
data indicate that this quantum critical point signals a 
quantum phase transition from the three-dimensionally 
ordered Neel phase to a spin-liquid state. Based on mag- 
netization data and an analysis of the microscopic details 
of the structure of BaCo2V20g, a collapse of the Neel or- 
der for this field direction has been predicted recently^ 
and is experimentally supported by our neutron scatter- 
ing data. 

Furthermore, our high-resolution measurements of 
the thermal expansion show that the Neel ordering in 
BaCo 2 V20g is accompanied by a structural transition 
from the tetragonal to an orthorhombic structure. This 
reduction of the structural symmetry can be straightfor- 
wardly explained by a finite magnetoelastic coupling of 
the inter-chain couplings and the fact that the magnetic 
symmetry is lower than tetragonal. The corresponding 
lattice distortions are rather weak and the macroscopic 
crystals are heavily twinned with two types of 90° do- 
mains. The domain orientation in the ordered phase can 
be influenced externally, either by applying uniaxial pres- 
sure parallel a t or by a magnetic field along the at axis. 
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